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Abstract. We construct an extension of the supersymmetric standard model where both CP symmetry and 
R-parity are spontaneously broken. We study the electroweak symmetry breaking sector of the model and 
find minima consistent with the experimental bounds on Higgs boson masses. Neutrino masses and mixing 
angles are generated through both seesaw and bilinear R-parity violation. We show that the hierarchical 
mass pattern is obtained, and mixings are consistent with measured values. Due to the spontaneous CP 
and R-parity violation, the neutrino sector is CP violating, and we calculate the corresponding phase. We 
further restrict the parameter space to agree with the limits on the electric dipole moment of the neutron. 
Finally, we study the CP violation parameter eK in the kaon system and show that we obtain results 
consistent with the experimental value. 

1 Introduction 

While the Standard Model of electroweak interactions (SM) has achieved a great deal of success, there still exist 
fundamental questions for which it does not provide answers. 

One of these challenging questions is charge-parity (CP) violation. Though at present the Cabibbo-Kobayashi- 
Maskawa (CKM) mechanism is in agreement with existing experimental data, it cannot accommodate the baryon 
asymmetry of the universe needed in big bang cosmology [1]. There is thus need for physics coming from beyond the 
SM scenarios. There are two basic possibilities to break CP: explicit at the Lagrangian level, or spontaneous in the 
vacuum [5]. The SM represents a case in which CP is broken explicitly [3 , through the introduction in the Lagrangian of 
complex Yukawa couplings which lead to CP violation in the charged-current weak interactions. Extensions of the SM 
which introduce new CP violating phases often lead to phenomenological difficulties. For instance, a general two Higgs 
doublet model with acceptable flavor changing interactions predicts a too large value for [?]. In supersymmetric 
(SUSY) extensions of the SM, one has additional sources of explicit CP violation arising from soft supersymmetry 
breaking terms in the Lagrangian. While in supersymmetric models a large number of new phases emerge [S], in a 
general minimal supersymmetric standard model (MSSM) these phases are strongly constrained by electric dipole 
moments (EDMs) [6]. Thus, in general, supersymmetric models share the problem of the origin of CP violation with 
the Standard Model and generate too large supersymmetric contributions to the dipole moments (known as the SUSY 
CP problem) |7j. 

An alternative scheme which could explain the source of CP violation is achieved through spontaneous symmetry 
breaking (SCPV) f2j. In this scenario the Lagrangian is invariant under the CP symmetry, but the ground state 
is asymmetric and the only sources of CP violation are the vacuum phases. Another motivation for SCPV arises 
from the strong CP problem. In spontaneous CP breaking vanishes at tree level and is calculable at higher orders 
|8]. Additional justification for studying SCPV comes from string theories, where CP exists as a good symmetry, 
but could be broken spontaneously [5]. Models with spontaneous CP violation require an extension of the minimal 
Higgs structure of the SM [10 . In general, more than one neutral Higgs particle will participate in flavor changing 
interactions. However, these interactions are severely restricted by the smallness of the Kj^ — Kg mass difference. One 
must have a mechanism to suppress such contributions either by extending neutral flavor conservation to the Higgs 
sector (requiring vanishing of flavor changing couplings) or requiring the neutral Higgs bosons to be heavy |llj . 
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Another problem of the Standard Model, which persists in the MSSM, is that the neutrino masses vanish. Yet the 
neutrino experiments have provided strong evidence for small nonvanishing neutrino masses |12j . Perhaps the most 
popular mechanism to explain neutrino masses is the seesaw mechanism [14,15 , which can generate small masses for 
neutrinos by allowing Majorana masses through the introduction of heavy right-handed neutrinos. Another popular 
way to explain the neutrino masses is through a small violation of R-parity [TS], Rp = (_i)3B-i+2s^ where i?=baryon 
number, L=lepton number, and s=spin of the particle |17j . Extensive phenomenological studies of R-parity violating 
effects exist; for a recent review see ^T. 

While there is no fundamental reason for the existence of R-parity in the MSSM, it is put in by hand in order to 
protect the proton from decaying. However models beyond the MSSM exist, which allow for R-parity violation in the 
lepton sector only. Such is the case if R-parity is broken spontaneously |18 l . If only lepton number (or baryon number) 
is violated, the proton does not decay. In, e.g., |19j . phenomenological implications of a model with spontaneous 
R-parity violation in the quark and leptonic sectors were explored. 

Both neutrino masses and CP violation could be explained if CP and R-parity were symmetries of the Lagrangian, 
but spontaneously broken by the vacuum. In this paper we consider a model where both of these violations are 
intertwined and both are spontaneous. CP is broken spontaneously through complex vacuum expectation values 
(VEVs) of the neutral scalar bosons, while R-parity violation and the seesaw mechanism will be allowed to provide 
the required masses and mixing parameters for neutrinos. We show that while it is non-trivial to satisfy conditions 
for both symmetries to break down at the same time, there are regions in the parameter space where suitable Higgs 
masses, as well as measured neutrino mass differences and mixing angles are obtained. In such a model, CP is violated 
in the neutrino sector. 

A model which included both neutrino mass generation mechanisms - namely seesaw and R-parity violation - 
with spontaneous R-parity violation was realized in .20J, where the spontaneous Rp violation was introduced via a 
term proportional to NLH2. This term represents the familiar bilinear Rp violating term mixing lepton and Higgs 
superfields, LH2, when the right-handed sneutrino field, N , develops a vacuum expectation value. This term, however, 
also breaks lepton number spontaneously, and thus introduces a superfluous massless Goldstone boson into the scalar 
spectrum. The problem can be solved by adding a singlet S to the theory, through a term N'^S, which explicitly breaks 
lepton number (but not R-parity), if S is not assigned lepton number —2. 

Attempts to violate CP spontaneously, by complex VEVs of the neutral scalars, exist ^ , but fulfilling the experi- 
mental constraints has proven difficult. More than one Higgs doublet is needed, see, e.g., [TTj. Spontaneous breaking of 
CP is not possible at tree level in the MSSM with two Higgs doublets, nor it is allowed in a model with four doublets 
[22] . Studies of minimal CP violations in the MSSM have shown that if no other symmetries are imposed, at least 
two extra singlet fields are required [3D]. Instead of adding doublets, or two extra singlets, one can study extended 
models, like the NMSSM model of [31j, where the so called /^-problem has been avoided by adding one singlet and 
requiring Z3 symmetry. At tree-level one cannot get spontaneous CP violation in this model either, and consequently 
radiative corrections must be evoked [32]. In that case a very light Higgs boson emerges [33] as it also happens in the 
MSSM, if spontaneous CP violation is induced via radiative corrections [33]. The consequences on the Higgs boson 
mass were also explored in |;35j- Another possibility studied is to discard the Z3 symmetry completely. On one hand, 
this way one loses the solution to the /i problem, but on the other hand, it is possible to achieve SCPV ^36j and also 
solve the problem of domain walls, which are created during the EW phase transition as the Z3 symmetry is broken 
spontaneously. 

An interesting model for spontaneous CP violation was presented in [37] , where the Z3 symmetry is replaced by R- 
symmetries on the whole superpotential, including non-renormalizable terms [38] . The method generates a Z3 breaking 
tadpole term for the singlet field S in the soft SUSY breaking part of the Lagrangian. This tadpole term allows for 
spontaneous CP violation to occur at tree- level [30II37| . The tadpole is assumed to originate from non-renormalizable 
interactions, which do not spoil quantum stability. We adopt this approach here. 

Several models in the literature have explored the consequences of breaking CP spontaneously. Some have studied 
the effects on leptonic observables, such as neutrino masses in the presence of right-handed neutrinos [23] and lepto- 
genesis [2T; as well as the effects on the electric dipole moments [25]. The consequences of having complex phases in 
the VEVs of the Higgs bosons have been analyzed in the kaon system [25] and B-meson system 27J. The consequence 
of allowing spontaneous, rather than explicit, CP breaking, is that the CKM matrix obtained is real. While several 
models mentioned above can generate CP violation in the kaon system that is consistent with the experimental data, 
a recent study argues that the CKM matrix is likely complex [2Hj; consequently a "hybrid" model was constructed 
in which more than one source of CP violation is present, allowing both a complex CKM matrix and non-trivial CP 
phases in the Higgs potential [2S]. Here we do not study the B-sector in detail and we expect that modifications in 
quark sector are needed, e.g., along the lines discussed above, to fulfill all experimental results. However, these changes 
do not qualitatively change the results obtained in this work, as we discuss later. 

Our paper is organized as follows. We give the Lagrangian and describe the model we used in Section [2] We explore 
Higgs boson masses and impose the condition that the masses satisfy experimental bounds in Section [3] We show that 
using both seesaw and R-parity violation, correct neutrino masses and mixings are obtained, in Section [4] We also 
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calculate the Jarlskog invariant of the neutrino sector. We explore the consequences of CP violation in the model and 
show that we can obtain a region of the parameter space compatible with the bounds on the electric dipole moments 
and obtain the observed ek in the kaon system in Section [s] We conclude in Section [6] 



2 The Model 



Our model is based on the superpotential 



2N, 



XhH^H^S) + '-^S' + '-^NfS, (1) 



where Hi and H2 denote the Higgs doublet superfields, Li and Qi the left-handed lepton and quark doublet superfields, 
respectively, and Ei and Ui, Di the lepton and quark singlet superfields. Right-handed neutrino superfields are denoted 
by Ni, i = 1,2,3, and S is the gauge singlet superfield. The SU{2) contraction is defined as £12 = —£21 = 1- Assuming 
baryon number conservation and Z3 symmetry, the terms in the superpotential are the only renormalizable ones that 
respect CP and R-parity, in addition to the gauge symmetry. All the parameters in the Lagrangian are real. 

The soft SUSY breaking terms in this model are the mass terms for scalars and gauginos, trilinear >l-terms, and 
the additional 5-tadpole, 



D*D, 



m2 



1 



MlS*S- 



{Msgg + M2WW + MiBB + h.c.) 



S'^ + £,^S + h.c. 



3! 



(2) 



Here i, j run over the family indices. In the tadpole term, £,^S, the parameter ^, which originates from nonrenormalizable 
terms [37], has been taken to be a free parameter, of the order of the soft supersymmetry breaking terms. We 
impose a flavor diagonal texture on h^j-,^ and the corresponding A-terms. The full tree-level scalar potential is 
Vs = Vsoft + Vf + Vd, where Vp and Vd are the usual F and D terms. All together, the model contains in the 
superpotential 51 additional parameters compared to the general MSSM. If we had instead the MSSM with additional 
explicit R-parity violating terms, there would be 48 new couplings T3j. If in addition CP would be broken, a large 
number of phases would appear in the soft masses and couplings. Compared to these numbers of parameters, our 
model with spontaneous CP and R-parity violation is economical. 

The minimization of the scalar potential with respect to the moduli of the scalar fields (jji and the corresponding 
phases 9i yield constraints later used in finding the scalar mass matrix. 



dVs 



0, 



dV 
d0,. 



(3) 



Without spontaneous CP violation, the VEVs are real and the minimization equations with respect to the phases 
are always satisfied. The minimization equations for the charged scalars can be trivially solved by setting all charged 
scalar VEVs to zero. As long as the tree-level masses of these fields remain positive and the corresponding soft A-terms 
remain small enough, this is also the global minimum of the potential with respect to these fields [2U] . The complex 
VEVs remain free parameters and we denote (the phase of Hi can always be rotated away): 



(Hi) 





V2e" 



(S) - ase" 



(4) 



Note that since Rp is violated, the W mass is to^ 



5(72^^ where 



al. w (174 GeVy 
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Minimizing with respect to the neutral scalars, we get at the tree- level: 

1^ = 2v,{m%. + X%{vl + al) + - vi + alj 

+ ta.nP{^XHXscrlcos{52 - 29s) + ^AffAjv.a^, cos((52 + + Ahcts 003(62 + Os))) 



2XHh%as(7L,(JR^ cos{9s - Ol, + Or,), 

+ h%ht{(TR^(TRk COs{0rj - Or^) + aLj(TL^ cos{Olj - Ol^)) 

+ cot (3{^XhXs(7s cos{52 - 20s) + ^Ai/AiViO-|. 003(^2 + 20^ 

+ 2A'^(TLi(JR. C0S(^2 + OLi - Or.) + 2h%XN^(TS<^L,(yRj COs((52 -Os + ^L^ + Or^), 



1^ = 2.2(ml,, + Xl{vl + al) - '-^{vf - vl + al) 
/(o-fl.o-Kfe cos{Orj -6r^)^ 
+ cot/J(^Ai/Asc7|cos(52 - 20s) + ^Ai/AiViO-|. 003(^2 + 20r,) + A^crs cos(52 + Os))) 



- — = 2as{ml + ^\{^\ + wf) + AsAhwiU2 cos((52 - 20s) + Aif^al 
das 

+ ]^Asasca3{2,0s) + ^A|a| + ^AsAjv,ct|. cos(26Is + 26iijJ) 

- 2^ cos{Os) + 2AhViV2 cos((52 + Os) + AN.aR, cos{Os - 20r,) 
+ 2h''^aLiaR. {Xrvi cos{Os - Ol^ + Or.) + Xn^V2 cos((52 -Os- Ol^ + Or.)), 

^ = 2a,AMl + ^^{v^ - vl + al)) 

+ 2A'^v2aR. cos((52 + OLi - Or^) + 2XHh'^viasaR- cos{Os - OLi + Or,) 
+ 2h^hi^vlaLj cos{OLi -Ol,) + 2h%XN^V2asaR^ cos(52 -Os + Or^ + 0^) 
+ 2h%h''rjaR^aL^aR, cos{OLi - Ol^ + Or, - Or^), 

dV 

= 2cri?,(M|_ + AN,ascos{9s - 29rJ + Xn,XhviV2 003(62 + 29rJ 

+ XN^Xsal cos(29s + 29rJ + X%al + XnM,'^% cos(29r^ - 9r.)) 
+ 2h^h^vlaR, cos(0R^ -Or,)+ 2A^>2CTL, cos(,52 - 9r^ + 9l^) 
+ 2XHh%viasaLj cos(9s + ORi -Ol^) + 2XN,h^NV2asaLj 003(62 - Os + ORi + Ol^) 
+ 2y^h^^aL,aR,^aL, cos(9r, - 9r, + 6l, - Ol,), 

— = XHViV2(-Xsas sin((52 - 29s) - ^Ni<y\ sm(62 + 26*^,)) - AHViV2as sin(52 + Os) 
002 

- 2v2(XNjh''^asaR.aLi sin(^2 - Os + Or. + Olj) + A^'^aR.aLi sin((52 - Or. + OLi)), 

dV 

—— = 2as(XHXsViV2as sin(^2 - 20s) - AhViV2 sm(6s + Os) + sm(Os) 
Cos 

- ^Asalsm(39s) ~ ^An^r^ sm(9s - 29rJ - ^XsX^^asa^ sm(29s + 29r^ 

- XEh'^viaR^aL, sm(9s + 9r. - 9l,) + XN^h'^V2aR^aL, sin((52 - 9s + Or^ + 9l,)), 
= 2aR,(aR,(AN.assm(0s - 20rJ - XhViV2 sia(62 + 20 rJ - lxsalsia(20s + 20 r,)) 

- (Ts(^Ljh%(XHVi sm(Os + Or. - OLj) + XNiV2 sin(52 -Os + Or^ + Olj)) 
+ A^NV2aL, sin(52 - Or, + Ol,) + h^h^vlaL, sm(OR>. - Or,) 

+ ^Ajv,Aiv,.CTfl,o-|. sin(26lfl. - 20r,) + h^^h'}}aR,aL^aL^ sm(eR, - Or. + Ol^ - Ol^)), 
—— = 2aLi(asaR h'^(XHVism(9s + Or. - OLi)- AAr,'(J2 sin((52 - Os + 9r + 6*^,)) 

OOLi 



- A%V2aR^ sm(62 - Or. + OLi) + ^jv ^^2^^!,^ sin.(OLk - OLi) 
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+ h%h)iaL^(TRMR, sm{9L, - 6*^. + Or, ~0r^)). (5) 

Here gi and 52 are the U(l) and SU(2) gauge couplings, respectively, and tan/3 = V2/V1. We use the seventeen 
minimization equations above solve for the soft masses of the neutral scalar fields and a subset of A-parameters 
{Ah, As, An,, A^^). 



3 Higgs Masses 

Separating the real and imaginary parts {(j) = (j^r + i4>i) of the nine neutral scalar fields (two Higgs, one singlet and six 
sneutrinos) we get an 18x18 dimensional mass matrix for the scalars. The radiative corrections to the scalar masses 
are implemented via the one-loop effective scalar potential [55] , 



Vl-loop 



327r2 



.5=1 ^ ^ a=l ^ 



(6) 



where m| are the field dependent eigenvalues of the 4 x 4 mass matrix and mf = — y^vf. 

renormalisation scale. The loop corrections lead to additional terms in both the minimisation conditions and the 
scalar mass matrix. In numerical calculations, we omit the D-term contributions and set for simplicity, Mq^^ = 
Mu^^ = MsusY, with MsusY ^ A^l TeV. 
The following experimental input is used: 

V ^ 174 GeV, mw ^ 80.42 GeV, mf "'"^ = 175 GeV, 

as = 0.102, mr = 1.777 GeV, sin^ 9^ = 0.23124. (7) 

Here 

Po/e 

™* = , *4 ■ (8) 

1+ 3^a. 

The rest of our free parameters are randomly sampled, with sampling ranges as follows (the couplings are constrained 
by perturbativity) : 

0.1 < Xh.n, < 0.4, 0.2 < As < 0.7, |/iAr| < 10"^ 
0.4 TeV < C < 1 TeV, -ir < 64, < tt, \{S)\ < 1 TeV, 

I (vl) I < 100 keV, I (N,) \ < 1 TeV, 2 < tan /3 < 60, (9) 



and the ^-parameters not eliminated by Eq. ([3| vary between < < (1 TeV)/i^. It should be noted that the 
VEVs of the right-handed sneutrinos are not constrained by any experimental bounds. In principle N could develop a 
VEV at a different scale than all the other scalars. However, if this VEV were to be complex, it would propagate CP 
violating phases to other (real valued) parameters of our model via their RG running equations and thus we restrict 
CT/j ~ MsuSY- 

In the limit Xn- ^ we recover lepton number conservation. CP and R-parity are still spontaneously violated as 
is lepton number, and thus in this limit the neutral scalar spectrum contains an additional Goldstone, as illustrated in 
Fig. [1] Using a model with only one right-handed neutrino (and consequently only one A at), the lightest scalar mass 
clearly tends to zero as the coupling Ajv ^ 0. 

For the full model we choose to set aR-^ ^ ~ Or^ ^ = to further reduce the sampling space. This choice identically 
solves the vacuum conditions dop^^ y = and de„^ V — Q. The ZZhi couplings {hi denotes any neutral Higgs) are 
reduced compared to the SM, due to other than Higgs doublet components in the physical Higgs bosons. Thus the 
Higgs strahlung production is less frequent than in the SM and experimental bound on the lightest Higgs mass is 
reduced from the SM value of ran > 114 GeV |40], see, e.g., discussion in [4r. Similarly the rate for the associated 
production of two Higgses through the couplings Zhihj is reduced. In Fig. l2] we plot the mass of the lightest Higgs 
boson as a function of tan f3. In the figure we have applied the experimental Imiits from LEP on all neutral spin-0 
particles to check that the masses are acceptable and we indicate the dominant component of the lightest Higgs. We 
have also studied the masses of the charged scalars at one-loop level and applied the experimental limits according to 
the main component of the charged scalar, i.e. if the main component is stau, we have applied the experimental limit 
of 81.9 GeV, as appropriate for stau |43]. Since the charged scalar can be relatively light, it is interesting to consider 
the possibility of seeing it at Tevatron in the decay t H^h. It appears that the light charged scalars are, however, 
mostly sleptons. Thus the coupling to quarks may be too weak to produce a significant branching ratio. 



6 



M. Frank, K. Huitu, T. Riippell: Higgs and neutrino sector in a spontaneous CP and R breaking model 



mjj(GeV) 



120 
100 
80 
60 
40 
20 



. 1 



0.2 

A 



0.3 



. 4 



N 



Fig. 1. Lightest scalar mass as a function of Aat. 
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Fig. 2. The lightest physical scalar mass as a function of tan (3. The main component of the lightest scalar is Hi 2 (orange 
rectangle), S (green diamond), ul (light blue triangle), or TV (dark blue cross). 



In Fig. |3] the composition of all the neutral Higgses is depicted as a function of their masses for one thousand 
parameter points satisfying the constraints mentioned above. It is seen that the light experimentally allowed Higgses 
tend to be mostly sneutrinos. In the region 100-150 GeV, a significant doublet Higgs component appears, as can also 
be seen from Fig. |2] showing that the lightest neutral Higgs is most often mostly doublet, if its mass is above 110 
GeV. The heavier Higgses are mostly either singlets or doublets. In Fig. |4]a Gaussian is fit to the number vs mass 
of the seven lightest neutral Higgses for one thousand parameter points. Curves for several of the lightest Higgses 
are strongly peaked, showing strong preference for particular mass values. The curves for heavier Higgses are much 
broader, showing much larger variation in their masses. Interestingly more than half a dozen are in the mass reach 
of the LHC. Unfortunately the doublet component in the light Higgses tends to be small, as seen in the Fig. ^ and 
their detection at LHC may be challenging. Detailed study of the detection is beyond the present work. 



4 Neutrino Masses and Mixing 



Since the scalars VEVs appearing in the neutrino mass matrix include phases, it is expected that the neutrino sector 
of the model is CP violating. 
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Fig. 3. Average composition of neutral Higgs particle as a function of the Higgs mass. Lightest grey/blue correspond to the 
sneutrinos, medium grey/green to singlet, and the darkest grey/orange to doublet Higgs. 




In a field basis of v^. ,Ni,S, Hi, B, W the neutral fermions form the following 11x11 mass matrix: 
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03X3 




03x1 


03x1 

















Olx3 




Xs{S) 








Olx3 










\h{S) 










Xh{S) 
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-72\^2 
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V2 



Ml 




v/2 \^L, 




M2 
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The mass matrix Eq. pO| differs from the neutrino mass matrix in [20] by the phases for the VEVs. It is easy to see 
the structure of the usual seesaw mechanism, which produces small neutrino masses m^, 







rriD 
Mr 



(11) 



where mu <C Mr. Similarly to |20j . there are actually several sources for neutrino masses: the usual seesaw and the 
mixing of neutrinos with and gauginos through R-parity breaking. 

It can be shown that the number of independent vectors in is an upper bound of the number of non-zero 
neutrino masses, e.g., in models with exclusively the gaugino seesaw, there is at most one nonzero neutrino mass at 
tree-level because there is only one linearly independent vector, in m^i. From (10 1 it is immediately apparent 



'Li I 



that we have four independent vectors in uiu '. three in the Yukawa matrix hjy , and the vector of sneutrino vevs 
If we were to include only one right-handed neutrino, there would be two linearly independent vectors, and as expected, 
we find that in such models there is one massless neutrino. In models with no right-handed neutrinos but bilinear 
R-parity violation, there are two independent vectors, /j.^ (= h]^ (N*) in our model) and The latter, however, 

can be rotated away using the accidental SU{A) symmetry of the {Li, Hi} fields, leaving only one independent vector 
and thus two massless neutrinos. 

Inspecting the requirement that rriD <^ A/fl'_yields some qualitative features of the model. In particular, the left- 
handed sneutrino VEVs must be small and h^iN*) should be of the same order. Thus, although {N*) is not bound by 
any other prior consideration, having /ijy ~ 10"^ results in an upper limit of a few TeV for the right-handed sneutrino 
VEVs. 

We diagonalise M^a numerically and use Mi ^ M2 ^ 1 TeV. Great care must be taken, as the elements of M^o may 
vary over ten orders of magnitude, and the eigenvalues themselves over as much as twenty orders of magnitude. Our 
calculations are carried out using a forced minimal precision of fifty decimals. The errors due to the lack of precision in 
this case appear farther than eight places behind the decimal point for the neutrino masses. The diagonalising matrix 



Af, with M* M^oAf- 



disLg^m^o , m^.) , has the following general form 



Af- 



(12) 



Here C,,C, <^ 1 denote 8x3 matrices that can be determined perturbatively, see e.g. [H]. Our interest lies in the 
matrix V^, the neutrino mixing matrix. Using the canonical notation for the neutrino mixing matrix [44], U = 
Vi, ■ diag(l, e'"^^, e*"^^) and c/s ij — cos / sin 9ij, 



C12C13 



S12C13 



sise 



Vy = I -S12C23 - Ci2S23Sl3e"^ C12C23 - 3125235136** S23C13 



S12S23 - Ci2S23Si3e 

we can extract the mixing angles as follows: 



i<5 



-C12S23 - Sl2C23Sl3e'* C23C13 



sin 6*13 = \ V^^ \ , tan 6*12 = 
One can also extract the CP violating Dirac phase 5 [15] : 

\5\— sin" 



V, 



11 



tan 023 = 



V, 



23 



cos ^13 sin 2^13 sin 2^12 sin 26*23 
where the Jarlskog invariant Jcp [46J of the neutrino sector is given by: 



(13) 



(14) 



(15) 



(16) 



In the quark sector, the J values in the SM are known to be J ^ 10 ^ . For points satisfying all the constraints from 
the scalar sector, we then apply the following experimental constraints concerning the neutrino sector |47j : 



sin^ 26*23 > 0.89, sin^ 6*13 < 0.047, sin^ 6*12 ~ 0.23 - 0.37, 



1.4xlO"^eV^ Q Q X. in-3„-i72 



7.3 X 10-^eV^ 



- 3.3 X 10"■'eV^ 
9.1 X lO^^eV^ 



(17) 



In Fig. |5] we show Jarlskog invariant in Eq. ( |16[ ) for a sample of 350 points in the parameter space satisfying all the 
scalar and neutrino sector constraints. It is seen that J is less than around 0.04. These values open the possibility of 
detecting the CP violation in the leptonic sector through neutrino oscillations, see [48ll49] and references therein. 
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Fig. 5. Jarlskog parameter as a function of the solar angle. 



5 Fermion Electric Dipole Moments and CP violation in the kaon system: ex 

Electric dipole moments represent a challenge for supersymmetric theories. It is known that the MSSM predicts too 
large EDMs by about three orders of magnitude for scalar fermion masses close to the current experimental bounds 
(0(100 GeV)) and CP violating phases of 0(1) [Tj- There are at present three solutions to this problem. One is to 
assume that supersymmetric phases are not of order unity, but rather of 0(10^^ — 10^'^) [7J. The second possibility is 
that the spectrum of the supersymmetric partners of quarks and leptons is heavy, i.e. of 0(3 TeV) or more |50j . and 
out of reach of the LHC. The third possibility is that there are internal cancellations among the different components 
of the neutron EDM (the chargino and gluino contributions in particular) which can reduce the magnitude of the 
neutron EDM [SP. Analyses have demonstrated that these cancellations are very difficult to achieve [ET . Finally 
attempts to set the flavor diagonal CP violation parameters to zero, but to allow CP violation through off diagonal 
elements in the scalar fermion mass matrices lead to too large EDMs, and further constraints must be imposed [53] . 
All these solutions are in effect fine tuning, either for the scalar fermion masses, or for the phases, or for part of the 
parameter space. 




Fig. 6. The loop contributions to fermion EDMs. 



The EDMs and in our model arise from loop contributions and are straightforward to calculate. The definition 
of the EDM for a spin-^ particle is 

Ci = ~'-df^a^,-f5WF'''', (18) 
and the general interaction Lagrangian between two fermions (W, ^) and a scalar (x) containing CP violation is 



ik 

This gives us the one loop EDM as 



1 - 75 , , 1 + 75 
k 7i r Lik 



if^Xfc + h.c.. 



E 



(47r)2m| 



QkB 



n2 \ n 



(19) 



(20) 
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Air) = ,^(3-r+'^), (21) 



2(1 - r)2 y 1 - r 

Bir)=:^y—r^{l + r+—^ . (22) 



2(1 -r)2 V 1-r 

Where Qi, Qk are the charges of the fermion and scalar respectively. Clearly we need lm{Kij^L*i^) for there to be 
a nonzero contribution. There are three different contributions depending on which particles are running in the loop: 
chargino, neutralino or gluino. As mentioned above, small EDMs are achieved if these contributions cancel out. From 



(20) the other two common solutions are also easily understood, since increasing the squark mass or suppressing 
im\KikL*i.) ^ both yield a small dj. 

The CP- violating parameter ex = ^'^{M.kk/'^'^k) receives no contribution from standard model processes in 
our model and the only contribution to M.kk comes from a chargino loop. Compared to the EDMs, the calculation 
of the kaon oscillation loop is more involved because there are non-perturbative hadronic states in the process. We 
use the vacuum insertion approximation (VIA) [llj whereby the matrix element {K\(srd){srd)\K) appearing in the 
loop calculation is reduced to two nonzero contributions (Vi and V2) which can be measured from kaon decays. 



— l- 





— 1 — ^ — 


4 1 




1 — 4 — 


S 1 




1 d 


s 




d 



Fig. 7. The loop contributions to eK- 



Defining the quark-squark-chargino interaction as 

- C,qqx = (Kj'fc7L + V^iffc7fl) qjXk + h.c. 
and introducing the following notation: 



Wl,u = Vj^,;vl^iV,%Vi, + {L^ R), 

(L^R), 
{L ^ i?), 



ijkl siL djL sjk aik 



Sti djL SJK dlK 

v^,*vI^,v^tvA 

^ siL djt * SJK * dik 



we can write the matrix element as 

1 

Mkk 



(27r)4 167r2 



-y 



ijkl 



_Tl 

2 ijkl 



-Vi - -V2 
3 2 
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(23) 



(24) 



Here I* are lengthy expressions arising from the loop integrals depending on the masses of the particles in the loop. 
Denoting Rab = {m^ — ml)/mf, they are 



Tl 

^ijkl 



T2 

'-ijkl 



jk 



R 



I ( {I — RikY 



ik 



log(l - R,k) 



l^{^~R^O? 



log(l-i?,,)- - 



^Ik 



i?-/(l - R,k) log(l - R,k) - Ra\l - R^l) log(l - Ra)] }• 



-^ik 



(25) 



(26) 
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Fig. 8. effect on neutron EDM. The bars show excluded ranges due to problems with vacuum stability (top), charged 

Higgs mass limits (middle), and neutrino sector constraints (bottom). 



Log I^K 




Ac (GeV) 



Fig. 9. eK as a function of the trilinear coupling A^. 



It is easily seen how these expressions simplify when using the mass insertion approxima tion (i.e. i = j). Since the 
model provides us with a full squark spectrum, we prefer to use the full expressions in Eq. (24). The VIA coefficients 
are given as follows [IT] : 



Vi 



K"''K 

2mK 



fK 160 MeV. 



(27) 



Using these results, our method to search for viable points in the parameter space has been to find first points 
which satisfy all the Higgs mass and neutrino sector constraints. From these points we begin a random walk in 
parameter space, slightly varying the parameter values for each step, checking the EDMs and discarding such steps 
as do not bring us closer to the experimentally acceptable values. There is a set of parameters {Au, M3, Ml/-) which 
enters exclusively in the calculation of the EDM. Suppression of the EDM can be achieved by increasing these mass 
parameters or, alternatively, one can find values for which cancellation between the different loop contributions occur. 
In general, varying only these parameters can lead to undesirably large mass parameter values (> 0(3 TeV). In such 
cases we vary all the parameters relevant to the EDM. Since the same parameters affect other sectors of our model 
as well, we have to check, at each random step, the various experimental constraints. In Fig. |8]we show the effect of 
changing Atv,, on the neutron EDM. The curve shows two values for Aat, where cancellations between the contributions 
to the EDM occur. It is seen that only one of the two dips in the curve satisfies all the required constraints. After this 
we try, again by randomly varying parameters, to find experimentally acceptable eK- This process is easier than for 
the EDMs, since the set of parameters {Ac, M'q , M^^), which only affects the value of ex, is sufficient for reaching 
acceptable ek values. In Fig. |9] the behaviour of ck, for the same parameter point as in Fig. [8] is shown as a function 
of one of the trilinear couplings, Ac. A clear peak where we find the experimentally allowed value can be seen. For the 
range shown, all the scalar and neutrino sector constraints, as well as EDMs, remain viable. 

We find that satisfying EDMs and kaon sector CP violation is possible without contstraining the phases, though 
quite restrictive, in the model. The reason why satisfactory values are found is due to the structure of the parameter 
space. There are variables, which affect EDMs or eK, but which do not affect the other observables. 
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6 Conclusions and discussion 

Violation of CP symmetry is well established, while neutrino masses make violation of R-parity attractive. We have 
considered here a model in which both R-parity and CP-symmetry break spontaneously. Our model contains in 
addition to the MSSM fields only the three right-handed neutrinos Ni and a singlet field S, which are both needed to 
spontaneously break CP and R-parity. We have shown that experimentally viable neutrino and Higgs boson masses 
can be obtained, and that CP is violated in the neutrino sector. We explored the parameter space requiring that EDM 
bounds and experimental results on the kaon system (especially ex) are satisfied. In [37] solutions for models satisfying 
both constraints in an R-parity conserving model were discussed. In our case, there are more parameters in relevant 
sectors, like the chargino mass matrix, which make the bounds mentioned easier to fulfill. 

While this paper is dedicated to the presentation of the model, and constraints on its parameter space, two 
immediate consequences arise for low energy phenomenology. In the leptonic sector, we predict a measurable CP 
violation, with the Jarlskog invariant below 0.04. In the scalar sector, we predict reduced ZZhi couplings (with hi 
any neutral Higgs boson) compared to the SM, and a Higgs strahlung production less frequent than in the SM. As 
well, the model favors a lightest Higgs mass less than its SM value, niH < 114 GeV. With our particle content, it is 
kinematically possible to produce several neutral Higgses within the LHC mass reach. 

In the present paper we have not studied the B-meson sector in detail, and we expect that modifications for 
the quark sector are needed before we can agree with the experimental results. Such modifications are achieved e.g. 
by using a model where the Higgs sector leads to a complex CKM matrix [23] or by adding vector quarks [29] . 
However, the parameters involved in the calculations of the neutrino sector and calculations involving quarks are 
largely disconnected, as we have seen from our kaon sector results. Thus the results of the present model concerning 
the Higgs and neutrino sectors are not expected to change qualitatively. The work along these directions is in progress 

This work is supported by the Academy of Finland (Project number 104368 and 115032) and by NSERC of Canada (0105354). 
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